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ABSTRACT 

A theory is developed which describes the local wave 
effects near a wave front advancing into calm water and is 
uniformly valid near the front. The waves generated by an 
oscillating vertical boundary (a wavemaker) and propagating 
out into an infinitely long channel are investigated both 
theoretically and experimentally. Comparison of the theo- 
retical and experimental results indicates that the theory 
does predict the wave phenomena about the front to a good 
degree of accuracy. 

Conclusions are drawn regarding the transient wave 
response to a wavemaker and compared to those of John W. 
Miles regarding the transient gravity wave response to an 
oscillating pressure. The transient features of the wave 
front considered are the width of the transient zone at a 
given time since the disturbance to the free surface of the 
fluid and the length of the transient period at a given 
distance from the disturbance. It is shown that the con- 
clusions regarding transient gravity wave phenomena apply 
without qualification to the waves generated by both oscil- 
lating pressure sources and moving surfaces. 
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NOMENCLATURE 



6- Dirac delta function 
n- Wave height 

g - Acceleration due to gravity 

G - Green's function 

T- Gamma function 

h - depth of fluid 

H - Heaviside step function 

I - Imaginary part of a complex variable 

R e - Real part of a complex variable 

s - Displacement of wavemaker paddle 

S - Amplitude of displacement of wavemaker 
paddle at free surface 

V - The ratio of distance to time (x/t) 

v - Same as V above. 

w - Frequency of oscillation of wavemaker paddle 
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INTRODUCTION 



The object of tnis study is to develop an expression 

wnicn describes the local wave effects near a wave front 

advancing into calm water that is uniformly valid about the 

front. Limiting forms of this expression as t -*■ 00 or x -*• 00 

were developed nistorically by Lamb'*' and by Cauchy and 
2 

Poisson respectively wnere t represents elapsed time after 
an initial disturbance to the free surface of the fluid and 
x the distance from the disturbance. Lamb's result de- 
scribes tne waves resulting from a steadily oscillating 
disturbance where the elapsed time since the disturbance 
was initiated is assumed to be so large that all transient 
features have disappeared and only the so-called steady 
state phenomena remains. The Caucny-Poisson expression 
describes the initial propagation of waves out from a dis- 
turbance but does not account for the time dependent aspects 
of tne disturbance. The exact formulation of these expres- 
sions depends on the nature of the disturbance, e.g., 
pressure variations or moving surface - and are widely 
available in the literature. 



'*'Lamb , H., "On Deep Nater Naves", Pvoc. Lond. Math. Soc., 
(2) , 2, 371, (1904) . 

2 

Cauchy, A., Poisson, S. D. "Memoire sur al theorie des 
ondes", Mem. de l' Acad. Roy. des Sciences , i. (1816). 



Neither of the above theories attempt to describe the 
transient wave phenomena spanning their respective limits. 
Tnis study develops a theory that predicts the waves in 
this transient zone and which upon passage to the appro- 
priate limit yields the nistorical result. The waves gen- 
erated by an oscillating vertical boundary (a wavemaker) 
and propagating out into an infinitely long channel are 
investigated both theoretically and experimentally. The 
channel is assumed to be infinitely wide so that the 
problem is two-dimensional. 

During the problem formulation phase of this analysis, 

tne attention of the author was directed to a paper by 
3 . 

Miles , which considered the transient gravity wave response 
to an oscillating pressure. Although the mechanism by which 
transient waves are generated is different in Miles ' 
approach, the qualitative features of the resulting waves 
snoulu be similar to tnose generated by a wavemaker. In nis 
paper Miles developed an expression predicting the transient 
waves tliac was uniformly valid about the wavefront. His 
solution method, though apparently consistent, was not 
immediately satisfying as it involved finding the asymptotic 
expansion of a double integral by applying the method of 
stationary phase twice in succession. It was felt that a 

> ----- a 

"’Miles, J. W., "Transient Gravity Nave Response to an Oscil- 
lating Pressure", j. Fluid Meo'h. t 13, (1962), pp. 145-150. 



more direct means of obtaining a uniformly valid expres- 
sion was desirable - both to verify Miles' result and to 
obtain clearer insight into tne transient problem. 

This study was undertaken primarily in the interest of 
academic completeness; viz., to bridge the gap between the 
solutions of Lamb and Cauchy-Poisson . However, several 
engineering applications can be foreseen. The utilization 
of wavemakers in towing tanks with corresponding applica- 
tions such as ship model and ocean structure testing is of 
direct interest to the naval architect or ocean engineer. 
Transient features of generated waves may allow some 
aspects of model testing to be explored more rigorously and 
completely. Further application is found in the area of 
ocean wave forecasting. The wave generating models used 
presently in predicting the wave phenomena associated with 
a storm at sea are not completely satisfactory and a new 
model, based partially or entirely on the expression 
derived in this study could be formulated. 



THEORETICAL ANALYSIS 



Consider the waves generated by a vertical plate in 
oscillatory motion given by (1) in a fluid bounded by a 
free surface, tne vertical plate, and of infinite extent 
otherwise. liiis so-called wavemaker problem is solved in 
Appendix B and the resulting expression for the waveheight 
as a function of distance and time is given by (2) . 



( 1 ) 



= S e W y/,g sin(wt) H(t) 



where S = amplitude of displacement at the 
free surface 

H(t) is the Heaviside step function 



( 2 ) 



n(x,t) = 



wS 



cos (ux) 
0 u + w 2 /g 



sin(wt) - sin(/gu t) 
w - /gu 



sin(wt) + sin(/gu t) 
w + /gu 



du 



Introducing the change in variable u = a 2 /g 



(3) n(x,t) = 



2wS 
" ) 



acos [ (a 2 /g) x] 
a 2 + w 2 



sin(wt) - sin(at) 



w - a 



sin(wt) + sin (at) 
w + a 



da = 



=1 



m 



2wS 



i (— x+at) i (— x-at) i (— x+wt) i ( u -x--wt) 



• ,o 



• ,a‘ 



• ,cr 
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a e ^ 



-^e g 



-wae 



+wae 



da 



Jo 



(a 2 +w 2 ) (a+w) (a-w) 



da 



a 



Using nondimensional variables x = w 2 x/g 
o/w and n = , (3) reduces to 



t = wt , 



(4) 



n(x,t) = 



r g 2 e i(g 2 x+qt)_ a 2 e i(g 2 x-at)_ ae i(a 2 x+t) +ge i(a 2 x-t) 
(a 2 +l) (g+1) (g-1) 



Equation (4) for n(x,t) may be considered as the sum 
of four integrals. 



(5) 



nU.U - i m {i x - i 2 - i 3 + x 4 } 



where 



x+at) 



v i (a : 
= f(g) e 

-'a 

= ff(a) 

J n 



da 



da 



1^ = 



I . = 



f(a) i(a 2 x+t) , 
— - — - e 'da 

a 



CUal e i(a 2 x-t) da 



and f(g) = 



g 4 -l 



The asymptotic behavior of these integrals for 
large x or t is readily evaluated by treating a as a com- 
plex variable and employing the residue theorem and 
Cauchy's integral theorem of complex variable theory. 
Their behavior for small x or t is given in Appendix H. 

In the complex o plane, the real integrals in (5) 
are transformed to complex integrals which are to be 



evaluated along the contour consisting of the ray from the 
origin to infinity along the positive real axis. 
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It is easily shown (see Appendix C) that the asymptotic 
behavior of integrals 1^, 1^ and 1^ is given by the follow- 
ing expressions for large t or x: 



^ if e i<x+t) 
I 3 -v, if e i(x+t > 

X il i (x-t) 

4 4 e 



+ 0 (t -3 ) 
+ 0(t~ 2 ) 
+ 0 (t~ 2 ) 



The asymptotic behavior of (5) then reduces to the 
following: 



( 6 ) 



n (x,t) 'v i m 



i Tt i (x-t) 

4 " 2 



The asymptotic behavior of is not as easily as- 
certained as was that of 1^, 1^ and I^. In the evaluation 
of these integrals, the best choice of an equivalent con- 
tour to facilitate the determination of their asymptotic 
behavior was obvious. If a similar contour is chosen for 
tne evaluation of I^, however, the integral diverges along 
portions of the contour. Consider then the exponential 

term of I 2 in expanded form where the substitution 
i 0 

a = r e has been made. It takes the form 
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i(o 2 x-at) it{r 2 (x/t) cos (29) -r sin(0)} x 
\ I ) © © 

x e ~t{r 2 (x/t) sin(20) -r sin(0)} 

A good choice for an equivalent contour would be one on 
which this term becomes small as x and t become large and 
vanishes when the contour is extended to infinity. For 
this term to become exponentially small for large values of 
t (or x) , the following restriction on r and 0 applies : 

(8) r 2 (x/t) sin(20) - r sin(9) > 0 

Let x/t = V: then for 0 negative r < 1/2V cos(0) 

for 0 positive r > 1/2V cos(0) 

for 0=0 r = 1/2V 

Thus , any contour with the general features of C 2 
shown below is a satisfactory equivalent contour for the 
path from 0 to + 00 along the real axis . The point 
a = 1/2V + Oi where the contour C 2 crosses the real axis is 
a stationary point of the integrand and must correspond 
to a saddle point as the integrand is a function of the 
complex variable 0 and therefore satisfies Laplace's 
equation . 
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a-plane 




Here singularities of tne integrand of have been 
indicated with x marks and in all cases are simple poles. 

The contour along the real axis is denoted by and has 
been indented at a = 1 . Thus, the above integral is 
defined in the Cauchy principle value sense. If the 
integral exists in the conventional sense, then its value 
will necessarily agree with the Cauchy principle value so 
there is no loss in generality. 

The existence of a saddle point at a = 1/2V suggests 
the most convenient equivalent contour on which to evaluate 
If C 0 is deformed so that it coincides with the 
steepest descents path from the saddle point, then divergent 



16 



integrals are avoided as the integrand of becomes 
exponentially small on such a path (granting suitable 
restrictions on f(o)). This contour is shown below. 
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Portions of C 2 that join the steepest descents path to 
the origin and to the positive real axis at infinity have 
been labeled and C R . The location of the saddle point 
relative to the simple pole at a = 1/2 V depends on the 
parameter V = x/t. Therefore, two general cases are shown. 

The closing contours and C R are chosen in the fol- 
lowing way: 

1) is the steepest descents path from the origin to 
infinity and intersects the contour at infinity. 

i 0 

2) C D is the arc of the circle Re joining the contour 
C 2 to the positive real axis. 

Consider the case V = x/t < 1/2. Application of 
Cauchy's integral theorem and the residue theorem lead to 
the following expression: 

(9) I 9 - i it Re s ( 1 ) = S r + S r + f r - 2iTiRes(l) - 2TriRes(-i) 

C 1 c 2 U R 

R-*-°o 



The asymptotic properties of each term in the above 
equality can now be evaluated. 



Res ( -i ) 



lim 
0-*— i 



(a+i)a 2 it(a 2 V-a) 
e 

a 4 -1 




-t -itV 
e e 



This term is exponentially small and can be neglected 
in comparison to terms algebraically small. 
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Res(l) = 



lim 

a-H-1 



(a-1) a 2 e it(a 2 V-a) 



a“-l 



]■ 



i(x-t) 



f r -* 0 as R -+■ 00 
R 



It is easily shown (see Appendix C) that f r = 0(t -3 ) 

C 1 

and is therefore small compared with other contributions to 
I 2 » Therefore, (9) reduces to 



(10) 



*2 ■ 



if e i(x - t > + 0 (t -3 ) 



Consider 



(11) 



it (c 2 V-a) 



C 2 c 4 -l 



da 



where V = x/t is a parameter to be specified and is less 

than 1/2 for the case under consideration. Introducing the 

1 2 

change in variable - ^ + ix = a V -a, the following con- 
siderations apply: 



4 ^ + ix = a 2 V - a 



y 2 



0 = 2V t 1 1 » / 4iV T 3 



2V 



% 



[ 1 + yt ] / Y = /4iV 



da 

dx 



+ Y 



4Vx 



Tz 
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On C 2 in the upper half plane and in the half 
plane to the right of C 2 



1 1 A 

0 = 2V [ 1 + YT 7 



YT /z ] = a + , 



da 



Y 



dx 4 Vt 



% 



On C 2 in the lower half plane and in the half plane 
to the left of C 2 

da_ ~y 



o = [ 1 - yi / 2 ]= o _ , 



dx 4Vx 



w 



Therefore, the integral along C 2 is given by 



( 12 ) 



a 2 it{a 2 (x/t) -a} 



c 2 a 4- i 



da 



- f -si- 

Jc„ a‘-l 



(a 2 V-a) 



da 



= e 



-it 

4V 



a 



2 da 



+ -tx , 

— e dx + 



o a + -l dx 



0 a 2 da 



- -tx , 

— e d x 



a 1 *-! dx 



If the Laurent expansion about x = 0 is obtained, the 
above integrals may be evaluated by integrating term by 
term. It is more convenient, however, to subtract off the 
pole at 0 = 1 and to then expand the remainder of the ex- 
pression about x = 0 so that the effect of the pole may be 
considered separately with a view towards obtaining a uni- 
formly valid result. These expansions are developed in 
Appendix D and are given below. 
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(13) 



a| da + -y/8 

o^-l dx 1-2V 




1-2V 






1-2V 



+ (-) {1-4V-4V 2 }t + 

(1+2V) (1+4V 2 ) 



y{ (1+4V 2 ) 2 -8V (1+2V) 2 } 
(1+2V) 2 ( 1+4V 2 ) 2 



+ (-) y 2 { (1+4V 2 ) 3 -4V(l+2V) 3 (3-4V 2 ) }t ^ , 
(1+2V) 3 (1+4V 2 ) 3 



(14) 



a ' 



da 



a -1 dx 



YA 



1-2V 



I, 1-2V “ 

x /z - 


-1 r 

Y 

+ — 


Y 


8 



-X 



-72 



1-2V 



{ 1-4V-4V 2 } 
(1+2V) (1+4V 2 ) 



t _1 /2 + y{ (1+4V 2 ) 2 -8V(l+2V) 2 } + 
(1+2V) 2 ( 1+4V 2 ) 2 



+ Y 2 { ( 1+4V 2 ) 3 -4V (1+2V) 3 ( 3-4V 2 ) } l /z + 
(1+2V) 3 (1+4V 2 ) 



Noting that terms containing integral powers of x 
cancel identically upon substitution into (12 ) , the 
resulting expression becomes 
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(15) 



g 2 e it(a 2 V-a) da _ -y/8 -it/4V 



2 0 4 - 1 




l- 


-2V 


/■°° “tX, 

e dx 


-it 






4V 

-4- Ck 




* 00 r 


ZVz 1-2 V 




4 


r\ V 


U l 

Y J 






U 



r 00 -tx , 
e dx 

V 2 . 12V 

0 X ' + 

Y 



y/8 x _ y/8 [1-4V-4V 2 ] 



1-2V (1+2V) (1+4V 2 ) 



y 3 /8 [ (1+4V 2 ) 3 -4V ( 1+2V) 3 ( 3-4V 2 ) ] x V 2 + Q ( T 3 / 2 
(1+2V) 3 ( 1+4V 2 ) 3 



-tx 



dx 



The first bracketed term in (15) includes the effect 
of the pole at a = 1. The second bracketed term can be 
expressed as a series of Gamma functions and is not in- 
fluenced by the pole at a = 1. Introducing the change of 
variable x /2 = s, the integrals in the first bracketed term 
above are transformed to the following: 



( 16 ) 



e ta dx 

J /?-l 1-2V 

Y 



,00 

o 



X 




dx 

1-2V 

Y 



,00 



-ts 



s + 



2sds 

1-2V 

Y 



+ 





2sds 

1-2V 

Y 



Let s = -s in the second integral of (16) . 



Then 
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(17) 



o° — ts 2 
e 2sds 

os + 1-2V 
Y 



-ts 



2sds 



r -ts 2 

' e 2sds 



s - 



1-2V 

Y 



s + 



1-2V 

Y 



0 -ts 2 



2sds 



s + 



1-2V 

Y 



J — o 



-ts ' 



s + 



2sds 

1-2V 

Y 



(s + i=™. - -ts 2 , 

^ y e ds 



s + 



1-2V 

Y 



= 2 



-ts 



ds - 



1-2V 

Y 



f 00 -ts , 
e ds 



J- c 



s + 



1-2V 
Y J 



y. 

With the change in variable y = t s the first inte- 
gral in (17) is readily evaluated and the second placed 
in a more convenient form. Equation (17) then reduces to 



(18) 



2 [ /w t _,/2 + f (±^> w(z) ] 



where W(z) = — 

ir 



e ^ dy 
z - y 



= - t V 2 ( lz2V) 

Y 



[I (z) > 0] 
m 



W(z) can be represented by the following expression : 



_ 2 

W(z) = e erfc(-iz) 



where erfc(z) = — 

/? 



e _t dt 



^Abramowitz, M. and Stegun, I. A. ed.. Handbook of Math, 
b'unj Lions j New York, (1965), p. 297. 
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On substituting for W(z), (18) can be placed in the 



form 



( 20 ) 



2 £ /? t 1/2 + l (i=~) W ( z) 



= 2 



, - ( 1 - 2V ) 2. 
/? t _/2 - iTi(i^X)e Y 



erf c U^)^}] 



Therefore, (11) can be expressed as follows for the 
case V = x/t < 1/2: 



( 21 ) 



c 2 it (a 2 V-a) t 
e do ^ 



a -1 



-it 



Zl/i ^ 4V 



1-2V 



, 1-2V v 2 

2iTi(i 2 ^)e Y erfc{i (— ^— ) t ^ 2 } 



„-it/4V | y/4 p ,l x ^-)/ 2 _ y/4 T ( 1/ 2 ) {l-4V-4V 2 )t ' 2 

1-2V 1 l 2 ; C " (1+2V) (1+4V) + 






+ zllZH d+4V 2 ) 3 ~4V(1+2V) 3 (3-4V 2 ) }T(3/2)t | 0(t - S / 2) 

(1+2V) 3 (1+4V 2 ) 3 



Substitution of (21) into (10) leads to the following 



expression for (6) : 
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( 22 ) 



n(x,t> * I m fee i,x - t) 



-it , 1-2V . 2 x- 

iTT ^ 4V 1 Y ' u ,1-2V. '/2 U 

— e e ' erfc{i( — - — )t } + 



+ e 



■ft 2 IT. 

- 1( 4? - ?> 



r V 



Vz 



{1-4V-4V 2 } _ i/ 

(1+2V) (1+4V2) r (V2)t 2 + 



+ ii t?y,) . 1 . ! . 3 . -4v 2 ) , ) r ( 3, t - v 2+ 0(t - 5 / 2) ll (v< i) 

(1+2V) 3 (1+4V 2 ) 3 A JJ Z 

where y = /4iV = V4i x/t 

Consider the case V = x/t > 1/2. Application of 
Cauchy's integral theorem and the residue theorem leads to 



(23) 



I_ - iiT Res(l) = + [ + [ 

z Jr Jr J ( 



<=L ; C 2 



R 

R-voo 



Similar considerations regarding 



and [ 

; C;l J ( 



R 



apply as for the case V = x/t < 1/2. Therefore, the 
asymptotic behavior of for large t or x is given by 



(24) 



I„ = 



it (a 2 V-o ) 



c 2 o“-l 



da + — e 



iir i (x-t ) , A . , -3 



+ 0 (t ) 



Proceeding as for the case V < 1/2, introduce the 
change in variable - ^ + ix = o 2 V - a. All previous con- 
clusions regarding the branches and branch cuts ofa(x) 
apply, hence 
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( 25 ) 



it (a 2 V-a ) 



da = 



C 2 a -1 



= e 



-it 

4V 



da 



— e *" T dx + 



da 



-tx , 

— e dx 



o a?-l dx 



o a -1 dx 



The Laurent expansions previously obtained about x = 0 
are valid, hence 




-it 






e 



4V 



y/4 r (1/2) t ^ 
1-2V 



y/4{l-4V-4V 2 }T (1/2) t 
( 1+2V) (1+4V 2 ) 



Y V 4 { ( 1 + 4V 2 ) 3 -4V ( 1+2V) 3 (3-4V 2 ) } r(3/2 )t' 3 / 2 + 0(t" 5/2 ) 
(1+2V) 3 (1+4V 2 ) 3 



-it 

Y/8 e 4V 
1-2V 




— ts 2 

e 2sds 
s + 1-2V 
Y 




-ts 



2 

2sds 



s 



1-2V 

Y 



The remaining integrals are treated in a similar manner 
as before except that the substitution s = -s should be made 
in the first integral in lieu of the second integral of (26) 
so that the resulting formulation exists. 
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(27) 



(28) 



-ts z , 
e sds 

“ 1-2V 
s + 



-ts ' 



sds 



o s - 



1-2V 

Y 





0 -ts 2 

e sds 


r°° -tS 2 , 

e sds 




J 


1-2V 


1-2V 


= 2 


_oo s - J 

V 


os- 

v 



— CO s — 



sds 

1-2V 

Y 



= 2 



-ts ' 



ds + 



1-2V 

Y 



■ts ' 



-°0 s - 



ds 

1-2V 



Let y 



= t*s. 



= 2 



/? 



t~ 1/2 -l (i^V) 



( ) W(z) 

i Y 



, , l /z , 1-2V) 

where z = t / ( — — — 



/? t >* + i7r( l^V )e ' Y 



_,1-2V. 2 

~ erfc{-i(^)t 1 /M 1 - 



J 



The evaluation of (25) is now complete and the wave- 



height n U, t) is therefore given by 



27 



(29) 



n(x,t) ^ { 



-it 2 / 1-2V . 2 

4r- e 4x e ^ er£c{ -i ( ^~ 2 ^ ) t } + 



+ e 



/t _ 2> i/ 

±( Jx 4 ; f V /2 /2{l-4V-4V 2 }r (1/2) - l / 2 



( 1+2V) ( 1+4V 2 ) 



+ 2iV^-l (1+4V 2 ) 3 -4V ( 1+2V) 3 (3-4V 2 ) } ,3. - 3 / 2 

(1+2V) 3 (1+4V 2 ) 3 2 



+ 







j 



(V > 1/2) 



where y = /4iV 

Tiie asymptotic: expansions for n(x,t) have been derived 
for the two separate cases of V = x/t < 1/2 and V = x/t>l/2. 
That these two expansions are equivalent can be shown by 
exploiting either of the following symmetry relations: 

(30) erfc(-z) = 2 - erfc(z) 

z 2 

(31) w(-z) = 2e - w(z) 

r+°° 2 

e /(z-t)dt (I z > 0) 



where w(z) = i/ir 



— oo 
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Considering equations (22) and (29) for r\(x,t), 
symmetry relation (30) is most convenient. 

From (30) , 



(32) 



erf c { i ( - ^ 2V -) t } = 2 - erf c {-i ( — - — ) t^ 2 } 



Substituting (32) into (22) one obtains 



(33) 



n (x, t) ^ i m • 



i tt i (x-t) 

e 

2 



in 

4 





1-2V 

Y 




erfc { -i ( 



1-2V ' 

Y 




+ 



+ e 



. ,t TT. 

1 4x ” 4 



(fV- 



(V < 1/2) 



Therefore 



(34) 



n(x,t) % i 



m 



-it 2 . 1-2V . 2 

” e 4x e * erfct-ili^X) t ^ ^ 



+ e 



• /t IK 

1 4x 4 



o (t 1/2 ) 



(V < 1/2) 



However, (34) is identical to (29) - the case V > 1/2. 
Therefore, the restriction V < 1/2 is unnecessary and can 
be removed and either expression completely describes the 
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waves. The final result then is uniformly valid about 
the wavefront whose position is defined by V = x/t = 1/2. 
Some simplification of n (x t) is possible and there are 
many equivalent expressions for it. A particularly compact 
expression is given below. 



ITT 



(35) 



n(x,t) ^ I m | ~ ^ erfc [ /x (l-t/2x) e 4 ] + 



+ 



. , t' TTv 

1 4x 4 



~ /i^7t{l-4x/t-4x 2 /t 2 }t 1/2 , 0(t - 3 / 2) 
2(l+2x/t) (l+4x 2 /t 2 ) 



For the cases V << 1/2, V >> 1/2, and V = 1/2, expres- 
sions for the waveheight ri(x,t) can be developed in an al- 
ternative manner. The effect of the simple pole at a = 1 
is either negligible or an easily evaluated residue contri- 
bution in these cases. These alternative expansions for 
n (x, t) have intrinsic merit and in addition provide a check 
on the validity of the uniformly valid expression that has 
been derived. They are available in Appendix E. 
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EXPERIMENTAL PROCEDURE 



Experiments were conducted in the Ship Model Towing 
Tank located in Building 48, Hydrodynamics Laboratory. The 
purpose of the experiments was to provide a comparison be- 
tween the waves generated by a wavemaker and those pre- 
dicted by the derived theoretical expression. 

Control of the oscillatory motions of the wavemaker 
paddle was accomplished by utilizing the sinusoid voltage 
output of a function generator as the input to the paddle 
motion control mechanism - a hydraulic device. A voltage 
output proportional to the deflection of the wavemaker 
paddle was fed back and compared to the input voltage pro- 
viding close synchronization of the paddle motion to the 
desired motion. The amplitude of horizontal displacement 
of the paddle at the water surface was measured. Both 
frequency of paddle oscillations and maximum amplitude of 
displacement were controllable parameters throughout the 
experiment. It was not found necessary to vary the maximum 
amplitude of displacement from its initial value of 1 inch, 
however. Paddle frequency was varied from 0.8 cycles per 
second (cps) to 2.0 cps . Above a frequency of 1.2 cps non- 
linear effects were encountered and severe vibrations of 
the wavemaker machinery noted. Transverse wave patterns 
were also induced indicating that significant three- 
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dimensional features were present in the flow. As a 
result, all experimental data for runs above 1.2 cps was 
not evaluated. Experiments were not conducted for fre- 
quencies below 0.8 cps so as to avoid the effects of 
finite towing tank depth. 

The height of the waves generated by the wavemaker 
was measured with waveprobes . These instruments are var- 
iable resistance devices whose electrical resistance changes 
with depth of immersion into the water. The voltage output 
of the waveprobes was connected to the input of a strip 
chart recorder (Sanborn, series 7700) . Thus a continuous 
record of the waves at a particular location from the wave- 
maker paddle versus time is obtained. Calibration of the 
waveprobes was accomplished by immersing the probes a 
measured distance into the water and adjusting the deflec- 
.tion of the strip chart recorder scribe to a convenient 
value. The deflection of the wavemaker boundary was also 
converted to a voltage signal and connected to the chart 
recorder input so that the motions of the wavemaker were 
continuously monitored. 
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RESULTS 



To facilitate comparison of the waves predicted by 
equation (35) and those measured experimentally, all dimen- 
sional quantities are normalized to a consistent base and 
the results illustrated graphically. The transformation 
equations relating dimensional to nondimensional quantities 
are those used in the theoretical analysis in the derivation 
of equation (35) and are as follows: 

x = w - — where w = frequency of oscillation 

g of wavemaker 

t = w t g = acceleration due to gravity 

— S = amplitude of oscillation of 

g = wavemaker at the free surface 

All bar quantities are dimensional 

Predicted values of the wave height are readily com- 
puted manually using tables to determine values of the error 
function complement for complex arguments. Equation (35) is 
well suited, however, for computer evaluation of n( x /t). 
Computer programs for this purpose are presented in Appendix 
G. In addition, a function subprogram for the error func- 
tion complement for complex arguments is listed separately. 

Computer facilities were also used to plot the theo- 
retical and experimental values of wave height using the 
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Stromberg-Carlson 4020 (SC-4020) cathode ray tube plotting 
device. Graphs of experimental results are computer faired 
with points corresponding to normalized, measured values 
of wave height encircled. Comparison of faired to experi- 
mental results show that the faired results differ very 
little from measured values. Reproductions of the wave 
records obtained during the experiments are available in 
Appendix F. 

Figures 1 through 5 compare theoretical and experimental 
results for five values of x. Wave height n(x,t) is plotted 
versus time t for fixed values of distance x. These figures 
therefore show the timewise distribution of wave height as 
seen by an observer located at a fixed distance from the 
wavemaker. Figure 6 is generated from equation (35) for 
x = 65 and illustrates well all the prominant features of the 
timewise distribution of waves. Figure 7 shows the spacewise 
distribution of wave height for various of time t. These 
figures are arranged sequentially and show the wave front 
propagating out from the wavemaker with increasing time. 



WAVEKEIGHT ■ WAV EH EIGHT - f|(x,t) 
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Fig. 1 -A . Wave Height vs. Time for x = 8.7687. Based 
on Equation (35). 




TIME - t 



Fig. 1 - B . Wave Height vs . • Titnejor x = 8.7687. 

on Measured Values (x = 11.6667 ft, 
f = 0.7826 cps). 



Based 



WAV EH EIGHT - " ( x , . ; WAVE HEIGHT - n(x,t) 
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Fig . 2-A . 



Wave Height vs. Time for x 
on Equation (35). 



13.6493. 



Based 



3.0 




Fig. 2-B . Wave Height vs. Time for x = 13.6493. Based 
on measured values (x ~ 11.6667 ft, 
f - 0.9764 cps). 



WAVEHEIGHT - ^(x,t) WAV EH EIGHT 
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3 .0 




Pig. 3-A. Wave Height vs. Time for x = 19.7632. Based 
on Equation (35) . 



a .n 




. J-B. Wave Height vs. Time_f or x = 19,7632. Based 
on measured values (x = 11.6667 ft, 
f = 1.1749 cps). 



Pig 



WAV EH EIGHT - ”>(x,t) WAV EHE I GHT - n ( x 
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Fig. 4-A. Wave Height vs. Time for x = 29.4239. Based 
on Equation (35). 




TIME - t 



Fig. 4-B. Wave Height vs. Time__for x = 29.4239. Based 
on measured values (x = 38.1667 ft, 
f = 0.7926 cps) . 



WAVEH EIGHT - n(x,t) WAVEHEIGHT 
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3 . 0 




Fig. 5-A. Wave Height vs. Time for x = 44.8634. Based 
on Equation (35), 



3.0 




Fig. 5-B. Wave Height vs. Time_f or x = 44.8634. Based 
on measured values (x = 38.1667 ft, 
f = 0.9787 cps) . 



WAVEHEIGHT - r|(x,t) 
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J.O 




TIME - t 

Fig. 6. Wave Height vs. Time for x = 65. Based on 
Equation (35) . 



Figure 6 illustrates well all the prominant features of 
the timewise distribution of waves . The wave height in- 
creases monotonically to a maximum value and then enters an 
oscillatory phase until the steady state wave amplitude is 
reached. 1 

Figure 7 shows the spacewise distribution of wave 
height for various values of time t. These figures are 
arranged sequentially and show the wave front propagating 
out into calm water with increasing time. 




WAVEHE IGHT - n(x,t) WAVEHEIGHT - r)(x,t) 
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5.0 




Fig. 7-A. Wave Height vs. Distance for t = 30.0. 
Based on Equation (35). 



5.0 




Fig. 7-B. Wave Height vs. Distance for t = 60.0, 
Based on Equation (35). 



WAVEHEIGHT - n(x,t) WAVEHEIGHT 



41 



S.O 




Fig. 7-C. Wave Height vs. Distance for t = 90.0 
Based on Equation (35). 




Fig. 7-D. Wave Heigh-t vs. Distance for t = 140.0. 
Based on Equation (35). 



WAVEHEIGHT - r|(x,t) WAVEHEIGHT - nCx.t) 
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Fig. 7-E . Wave Height vs. Distance for t => 180.0. 
Based on Equation (35). 




Fig . 



Wave Height vs. Distance for t = 
Based on Equation (35). 



7-F 



240.0 



I 














DISCUSSION 
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The waves predicted by the equation (35) for the wave 
height as a function of distance and time compare well with 
measured values. To an observer located a fixed distance 
from the wavemaker, the wave height increases monotonically 
with time to a maximum value and then enters an oscillatory 
phase during which the steady state wave amplitude is 
approached. Both measured and predicted wave records show 
this characteristic behavior. 

Measured and predicted values of the wave height cor- 
respond very closely in phase, except for very small values 
of t when equation (35) is no longer valid. However, meas- 
ured values of wave amplitude are consistently larger than 
predicted values. This difference can be attributed to the 
boundary condition imposed on the moving surface of the 
wavemaker in the theoretical development. For convenience 
of solution, this boundary condition was assumed to be the 
following : 

(1) s = S e W ^sin(wt) H(t) (-°° < y < 0) 



However, a more realistic approximation to the flow 
near the towing tank wavemaker paddle is given by 



s 



S ( ? + ^ sin(wt) H(t) 






(y _> -h) 



(36) 



0 



(y < -h) 
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This is not the exact boundary condition but is a reason- 
able approximation to it. It assumes that displacements of 
the wavemaker from the vertical are small with respect to 
the depth of the wavemaker paddle (approximately equal to 
the depth h of the towing tank wave generator) and is the 
linear approximation to the actual displacement. 




In addition, the theoretical expression for ri(x,t) is for 
the case of infinite depth water. The waves generated by 
the towing tank wavemaker have been chosen so that finite 
depth effects have a negligible effect on the wave features. 
Therefore, the boundary condition given by (36) is a simpli- 
fied model of the physical situation shown above. 

A theoretical development based on the more realistic 
approximation to the boundary condition poses no additional 
complication to the problem but does significantly add to 
the algebraic tedium and perhaps obscures evaluation of the 
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results somewhat. Of primary concern is the determination 
as to just how the use of boundary condition (1) above in 
lieu of (36) has influenced the problem. Physical consid- 
erations suggest that assuming an exponential decrease in 
fluid velocity with depth instead of the more realistic 
linear decay should result in wave amplitude predictions of 
smaller magnitude than measured values. As less energy is 
imparted to the fluid near the moving surface (average 
amplitude of oscillation being smaller) , wave amplitudes 
must be correspondingly smaller as energy considerations are 
ultimately related to the square of wave amplitude. This is 
substantiated by experimental results. 

Partial verification of the physical argument presented 
above can be readily obtained in spite of the algebraic 
complexities of the more realistic boundary condition. The 
steady state wave amplitude can be easily obtained by let- 
ting t approach infinity and evaluating the remaining terms. 
For the case at hand this amounts to determining the residue 
contributions to the integral representation for the wave 
height . 

The solution to the unsteady wavemaker problem (in- 
finite depth water) can be placed in a more general form 
where the functional dependence of the wavemaker paddle 
displacement on depth is left unspecified (subject to the 



restrictions previously mentioned on this function) . One 
formulation of it is the following: 
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(37) 



n(x,t) = - ^ - 



cos (ux) 



e Uy F (y ) dy 



sin (wt) -sin ( /gu t) 
w - /gu 



+ sin(wt)+sin(/gu t) 1 du 
w + /gu -* 



For F(y) = 



w 2 

T y 

e s , 



e uy F (y) = - 



u + w 2 /g 



and equation (2) is rederived, 



y + h 
h 



For F(y) = < 



y >_ -h 
y < -h 



e Uy F(y)dy = - — ■ (uh + e - 1) 
0 hu 2 



Therefore, the wave height n (x,t) is given by 



(38) n (x,t) = 



2ws 



-uh 



(uh+e 
hu 2 



1) , , T/gu sin(/gu t)-w sin (wt)~l , 

— cos (ux) — 3 — 3 - du 

gu - w 2 
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Introducing the change in variable u = a 2 / g and 
using non-dimensional variables x = w 2 /g x, t = wt , a = o/vj, 
h = w 2 /g h, and r| = Trri/2s one obtains: 



(39) 



n (x, t) = i 



m 



{a 2 h+e a ^-1} / „ i(a 2 



ha 4 (a+1) (a-1) 



ae 



! x+at) 



ae i(a 2 x-at) _ e i(a 2 x+t) + e 



i(a 2 x-t) | 



da 



Following methods similar to those used in deriving 
(35) , as t + “ 

(40) n (x,t) 'v- I m { iir Res[l] } 



a. 



(h + e ^ 
h 



1) TT 



cos (x-t) 



where h = depth of towing tank. 

For the measured depth of 43.5 inches (3.62 feet) and 
wavemaker paddle frequency of 0.9787 cps , h = 4.25. 



(41) n(X/t) ^ 



(4.25 + e 4,25 -l 
4.25 



X 7T COS (X-t) 



2 . 4 cos (x-t) 



Equation (41) compares favorably with measured steady 
state values of wave amplitude corresponding to a paddle 



frequency of 0.9787 cps . These experimentally obtained 
values of wave height are given by: 
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(42) 



ri(x,t) - 2.35 cos (x-t) 



The ratio of the steady state wave height associated 
with boundary condition (36) to that associated with 
boundary condition (1) is 2.4/1.57 = 1.53. If this ratio 
is assumed to hold in a very loose sense throughout the 
transient zone, then the predicted values of wave height 
are very close to measured values. 

For large values of t or x, equation (35) for r|(x,t) 
must approach the historical solutions of Lamb and of 
Cauchy and Poisson. This property can be shown by expanding 
the erfc for large positive and negative values of its ar- 
gument. Letting t approach 00 one obtains: 

(43) r|(x,t) 'v ir/2 cos (x-t) 

These are the outgoing waves at infinity that Lamb derived. 
Letting x approach <» one obtains as the first order term: 



(44) 




Equation (44) exhibits the general features of the 



Cauchy/Poisson solution. 
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CONCLUSIONS 



In his analysis of the transient gravity wave response 



qualitative features of the generated waves as were observed 
in the towing tank and predicted by the derived expression 
for n (x, t) , Equation (35). The wave amplitudes associated 
with a concentrated oscillating pressure also show the mono- 
tonic increase to a maximum value followed by an oscillatory 
phase to the steady state amplitude. 

In addition. Miles defines a rise-time T as the time 
for the timewise envelope to rise from 0.1 of its steady 
state value to its first maximum, and a width X of the 
transition zone for the spacewise envelope as the distance 
between the point where the wave height reaches 0.1 of its 
steady state value and its maximum amplitude. The expres- 
sions for T and X are as follows : 



It is of interest to determine if the theory developed for 
the transient wave response to a wavemaker supports these 
observations or modifies them in any meaningful way. 



to an oscillating pressure. Miles'* described many of the 



(45) 



T - 9 (x/g) ‘ 



(46) 




5 



Miles, Op. Cit. , pp. 149-150. 
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Equation (35) gives the transient wave response to a 
wavemaker and can be rewritten as follows : 

(47) n(x,t) 'v I m e 1 ^ erfc [/x (l-t/2x) e i77//4 ] + 0 (t ^ 2 ) 

Terms of order 0(t /2 ) can be neglected for sufficiently 
large t (or x) . It is convenient to treat the remainder of 
(47) as an outgoing plane wave with amplitude and phase 
given functions of x and t. With this view in mind, (47) 
can be rewritten as follows: 



(48) 


n (x. 


t) 'V A cos [ (x-t) + ifl] 


where 






(49) 


A = tt/4 


| erf c [ /x ( l-t/2x) e"' -77 '^ 4 ] | 


(50) 


ip = Arg 


{ erf c [ /x (l-t/2x) e"*" 77 ^ 4 ] } 



Expressions corresponding to (45) and (46) can be 
developed for the transient wave response to a wavemaker 
by investigating the behavior of A and using Miles 1 
definitions of X and T. 

The behavior of A as a function of U =/x(l-t/2x) is 



shown below. 
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5 




Fig. 8. Timewise Wave Envelope 



A is equal to 0.1 of its steady state value (A = 0.15708 = 
0 . lx n/2) when U = +2.79 5. A reaches its maximum value when 
U = -1.5255. With these values, the rise-time T and the 
width of the transient zone X are readily determined. In 
terms of x and t, the maximum and 0.1 steady state ampli 
tudes occur when: 



(51) 


/x (l-t/2x) = -1.5255 


for A = A 

max 


(52) 


/x (l-t/2x) = 2. 795 


for A = 0.1 xtt/ 2 = A Q>lss 
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(53) 



For x = x : 
o 



t m = 2 x q + (2) (1.5255) x q ' 2 

t = 2 x q - (2) ( 2 . 7 9 5 ) x ^ 2 

O.lss ° 



Therefore 



T = t - t n = 8 . 64 lx 

m O.lss o 



Vz 



T - 9x 



y 2 



For t = t : 
o 



/x (1-t /2x ) 
m o ra 



= -1.5255 for A = A 



max 



^0. 1 ss (1 -V 2x 0.1ss ) - 2 ' 755 for A = fl o.lss 



These relations for x reduce to 

m 



x 2 - (t +(-1.5255) 2 }x + t 2 /4 = 0 

m o m o 



x 0 . lss - (t o +(2 - 795)2)x 0.1ss + fc o /4 - 0 

x and x_ . are then given by 
m O.lss J J 

{t Q + (-1. 5255) 2 } - 1.5255/2t +(-1.5255) 2 



{t +(2.795) 2 } + 2. 795/2t +(2.795) 2 
_ o o 

x 0 . lss 



2 



